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Abstract. We study linear integro-differential equations in Hilbert spaces with 
operator-valued kernels and give sufficient conditions for the well-posedness. We 
show that several types of integro-differential equations are covered by the class 
of evolutionary equations introduced in [R. Picard. A structural observation for 
linear material laws in classical mathematical physics. Math. Methods Appl. Sci., 
32(14): 1768-1803, 2009]. We therefore give criteria for the well-posedness within 
this framework. As an example we apply our results to the equations of visco- 
elasticity. 
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1 Introduction 



It appears that classical phenomena in mathematical physics, like heat conduction, wave prop- 
agation or elasticity, show some memory effects (see e.g. |1HI16|). One way to mathematically 
model these effects is to use integro-differential equations. In this work we give a unified ap- 
proach to the well-posedness of linear integro-differential equations of hyperbolic and parabolic 
type, i.e. equations of the form 

t t 
u{t)+ j g{t - s)u{s) ds + A* Au{t) - j A* h{t - s) Au{s) ds = f {t) (t G M), (1) 

— oo — oo 

and of the form 

t t 

u{t)+ j g{t - s)u{s) ds + A* Au{t) - j A*h{t - s)Au{s) ds = f{t) (t G M), (2) 

— oo — oo 

respectively. In both cases A denotes a closed, densely defined linear operator on some Hilbert 
space, which is in applications a differential operator with respect to the spatial variables. 
These type of equations were treated by several authors, mostly assuming that the kernels 
g and h are scalar-valued, while we allow g and h to be operator-valued. If the kernels are 
absolutely continuous, the well-posedness can easily be shown. However, we focus on kernels, 
which are just integrable in some sense without assuming any kind of differentiability. 
The theory of integro-differential equations has a long history and there exists a large amount 
of works by several authors and we just mention the monographs |2H I12| and the references 
therein for possible approaches. Topics like well-posedness and the asymptotic behaviour of 
solutions were studied by several authors, even for semi-linear versions of ([1]) or ([2]) (e.g. 
O [21 H] for the hyperbolic and jl] El [3] for the parabolic case). 

Our approach to deal with integro-differential equations invokes the framework of evolutionary 
equations, introduced by Picard in [18| 119) . The main idea is to rewrite the equations as 
problems of the form 

{doM{d^^) + A)U = F. (3) 

Here do denotes the time-derivative established as a boundedly invertible operator in a suitable 
exponentially weighted L2-space. The operator M(5(^^), called the linear material law, is a 
bounded operator in time and space and is defined as an analytic, operator-valued function of 
. The operator A is assumed to be skew-selfadjoint (this can be relaxed to the assumption 
that A is maximal monotone, see I26|). As it was already mentioned in |15| . the operator 
M (Sf)"^) can be a convolution with an operator-valued function and we will point out, which 
kind of linear material laws yield integro-differential equations. By the solution theory for 
equations of the form ([3|) (see |18| Solution Theory] or Theorem 12.91 in this article) it suffices 
to show the strict positive definiteness of d^M^dQ^) in order to obtain well-posedness of the 
problem. Besides existence, uniqueness and continuous dependence we obtain the causality of 
the respective solution operators, which enables us to treat initial value problems. 
In Section 2 we recall the notion of linear material laws, evolutionary equations and we state 
the solution theory for this class of differential equations. Section 3 is devoted to the well- 
posedness of hyperbolic- and parabolic-type integro-differential equations. We will show how 
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2 Evolutionary Equations 



to reformulate the problem as an evolutionary equation and state conditions for the involved 
kernels, which imply the positive definiteness of y{tdoM{dQ^) and therefore yield the well- 
posedness of the problems. Furthermore, in Subsection 3.1 we will briefly discuss a way 
how to treat initial value problems (see Remark I3.10p as well as problems where the whole 
history of the unknown is given (Remark 13. lip . Finally, we apply our findings in Section 4 
to the equations of visco-elasticity. This problem was also treated by Dafermos [HI [7], even 
for operator-valued kernels but under the stronger assumption that the kernels are absolutely 
continuous. 

Throughout, every Hilbert space is assumed to be complex and the inner product, denoted by 
{■\-)h is linear in the second and anti-linear in the first argument. Norms are usually denoted 
by I • I except the operator- norm, which we denote by || • ||. 



In this section we recall the notion of evolutionary equations due to |18| [T9l [20] . We begin 
to introduce the exponentially weighted L2-space and the time-derivative do, established as 
a normal, boundedly invertible operator on this space. Using the spectral representation 
of this time-derivative operator, we define so called linear material laws as operator-valued 
functions of d^^ . In the second subsection we recall the solution theory for evolutionary 
equations [18| Solution Theory] and the notion of causality. 

2.1 The Time-derivative and Linear Material Laws 

Throughout let z/ G M. As in |20[ \TE\ [T5] we begin to introduce the exponentially weighted 
L2-spac^. 

Definition 2.1. We define the Hilbert space 



2 Evolutionary Equations 



H^fi(^) := if :R^C f measurable, / |/(t)|e 




endowed with the inner-product 




Remark 2.2. Obviously the operator 



e 



: i7,,o(M) ^ L2{R) 



defined by {e-"""/) (t) = e'^^fit) for t G M is unitary. 



^For convenience we always identify tlie equivalence classes with respect to the equality almost everywhere 
with their respective representers. 
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2.1 The Time-derivative and Linear Material Laws 
We define the derivative d on -L2(M) as the closure of the operator 

(j) ^ (/)', 

where (M) denotes the space of infinitely differentiable functions on M with compact sup- 
port. This operator is known to be skew-selfadjoint (see |29[ p. 198, Example 3]) and its 
spectral representation is given by the Fourier- Transform J^, which is given by 

{Tcj)) (t) := [ e-'"V(s) ds {t G M) 



for functions (p G Li(M) n -L2(M), i.e., we have 

d = T*{im)T, (4) 

where m : D{m) C L2(K) — )• L2(ffi) denotes the multiplication-by-the-argument operator 
{{mf) (t) = tf(t)) with maximal domain D(m). 

Definition 2.3. We define the operator on H^fi{M.) by 

and obtain again a skew-selfadjoint operator. From ^ we immediately get 

which yields the spectral representation for dy by the so-called Fourier- Laplace- Transform 
U ■■= J'e-^"' : Hy,o{R) L2(M). 

An easy computation shows, that for (p G C^(M) we get cp' = d,y4> -\- which leads to the 
following definition. 

Definition 2.4. We define the operator dQ^y := dy-\-u, the time- derivative on Hyfi{K). If the 
choice of G M is clear from the context we will write instead of d^^y. 

Remark 2.5. Another way to introduce 9o,v is by taking the closure of the usual derivative of 
test-functions with respect to the topology in Hyfl{W)., i.e. 

-^^■.,o(K)e^^'.,o(IR) 



We state some properties of the derivative (9o,i/ and refer to |15| [20] for the proofs. 
Proposition 2.6. Let u > Then the following statements hold: 

(a) The operator do^y is normal and G Q{do^y) with WOqIW = K 

(b) do^y = C*y(\m + v)Cy. 
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2 Evolutionary Equations 



(c) For u € HyfiiM) we have (t) = u{s) ds for almost every t G 

Of course, the operator d^^i, can be lifted in the canonical way to Hilbert space- valued functions 
and for convenience we will use the same notation for the derivative on scalar-valued and on 
Hilbert space- valued functions. The space of Hilbert space- valued functions, which are square- 
integrable with respect to the exponentially weighted Lebesgue measure will be denoted by 
Hpfl{M.] H) for 1/ S M. Using the spectral representation for the inverse time-derivative d^l, for 
> 0, we introduce linear material laws as follows. 

Definition 2.7. For r > let M : Bc{r,r) — >• L{H) be a bounded, analytic function. Then 
we define the bounded linear operator 

M (^^—] : L2{R;H) ^ L2{R;H) 
\im + v J 

for 1/ > ^ by i^M (^j;^^ {t) = M (^jj^^ f{t) for t G M and the linear material law 
Mid^^l) by 

M{do^l) := CIM (-^) G L{H,^o{R;H),H,^o{^;H)). 

Note that the operator Mip^]^^ as a function of d^\,, commutes with the derivative d^^y^ in 
the sense that d^^yMidQ^) D M{dQ l)do,u- 

Remark 2.8. The assumed analyticity of the mapping M is needed to ensure the causality (see 
Theorem 12. 9p of the operator M((9(^jJ,) using a Paley- Wiener- type result (cf. |23)). 



2.2 Well-posedness and Causality of Evolutionary Equations 

In |18) the following type of a differential equation was considered: 

(^do,uM{d^l) + a) n = /, (5) 

where A : D{A) C H ^ H is a, skew-selfadjoint operator, / € Hiy^Q{M;H) is an arbitrary 
source term and u G Hi^^q{R; H) is the unknown. For this class of problems the following 
solution theory was established. 

Theorem 2.9 ( |18[ Solution Theory]). Let A : D{A) <^ H ^ H be a skew-selfadjoint operator 
and let M : Bc{r,r) — ?• L{H) be analytic, bounded and such that there exists c > such that 
for all z E B£(r,r) the following holds 

mcz-^Miz)>c. (6) 

Then for each ^ > the operator ^5o,i/M(9(^jJ,) + A^ is boundedly invertible as an operator 
on H,^fi(E.; H) and the inverse is causal, i.e. 

XR<S^){do,,M{do^l) + A) = XR<„(m)(5o,.M(5o-^) + ^)"' 

■^This shows, that for u > the operator d^l is causal, while for ly < we get the anti-causal operator given 
by Oq Iu = — J°° u{s) ds (see [15]). 
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for each a G mJI 



This means that under the hypotheses of Theorem I2.9[ Problem ([5]) is weh-posed, i.e. the 
uniqueness, existence and continuous dependence on the data / of a solution u is guaranteed. 
However, ([5]) just holds in the sense of 



do,uM{d^l)+A)u = f. 



where the closure of the operator is taken with respect to the topology on Hiyfi{M.; H). To 
avoid the closure, one can use the concept of extrapolation spaces, so-called Sobolev-chains 
with respect to the operator A + 1 and (9o,jy (see |17| . |20| Chapter 2]). In this context 
Equation ([5]) holds in the space H,^^^i{M; H-i{A + 1)), where we denote by (i7jy^fc(R))^g2 the 
Sobolev-chain associated to (9o,i/- Using that M{dQl) and A commute with Oq,^, one derives 
the following corollary from Theorem | 



Corollary 2.10. Under the conditions of Theorem \2.y\ the solution operator \^dQ^i^M{dQ I) + A 
extends to a bounded linear operator on //^. ^(IR;//) for each k £ 

Remark 2.11 ( |20| Chapter 6], [281 Theorem 1.4.2]). The solution theory is independent of the 
particular choice of > ^ in the sense that for right-hand sides / S H,y^k(^', H) n //^^^(M; H) 
for n^v > A; G Z we have 

d,,,M{d^l) + a)"' / = (5o,mM(5o-;) + a)~^ f. 



3 Integro-Differential Equations 

In this section we introduce an abstract type of integro-differential equations with operator- 
valued kernels, which covers hyperbolic- and parabolic-type equations. This abstract type 
allows to treat convolutions with the unknown as well as with the derivatives (with respect to 
time and space) of the unknown. We introduce the space Li^^{M>q; L{H)) for fi £ M as the 
space of weakly measurable functions B : M>o — )• L{H) (i.e. for every x,y £ H the function 
1 1— >■ {B{t)x\y) is measurable) such that the function 1 1— )• ||S(t)|| is measurabl^ and 

oo 

I^Il,„(m>o;L(//)) := I e-^'\\B{t)\\ dt < oo. 



Note that Li,^(M>o; LiH)) ^ Li,,,(M>o; L{H)) for fi < u. For a function B G Li,^(M>o; L{H)) 
we can establish the Fourier-transform of i? as a function on the lower half-plane [M] — i[M>^] := 
{t — ii^ 1 1 G R, > //} by defining 

oo 

{B{t - ii^)x\y) := [ e-'^'e-"' {B{s)x\y) ds (t G M, > /z) 

v2vr J 


^Here we denote by XK<a('Ti) the cut-off operator given by (xK<a('^)/) (*) ~ Xa<„ 

*!£ H is separable, then the weak measurabihty imphes the measurabhtiy of t i— >■ ||_B(t)||. 
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3 Integro-Differential Equations 



for x,y £ H. Obviously the function t — iz/ i— )• B{t — iu) is bounded on [M] — i[M>^] with values 
in the bounded operators on H and satisfies |-B|i^([K]_i[R>^,];L(H)) < ^l^lLi,f,{R>o;L(/f))- 
Moreover it is analytic on the open half plane [M] — i[R>^]2|. For B G Li^^{M.>q; L[H)) we 
define the convolution operator as follows. 

Lemma 3.1. Let B G Li^^(M>o; for some ^ G M. We denote by S{W;H) the space of 
simple functions on R with values in H . Then for each f > fi the convolution operatoi0 

B* : 5(M; H) C H^^oi^; H) i/^,o(M; H) 



u\-^ \t^ I B{t - s)u{s) ds 
\ I ) 



is hounded with \\B * \\l{h^o{R;H)) — „{R>o;L(H))- Hence, it can be extended to a bounded 

linear operator on H^fi(M; H). 

Proof. Let u > ^. Then we estimate for u G S'(M; H) using Young's inequality 

e-2'^*|(5* n)(t)|2 dt < / I / e-^(*-")||S(t-s)||e"'^^|n(s)| ds I dt 



< I J e-'''\\B{t)\\ dt\ J |n(t)|2e-2'^* dt, 



which yields B * u £ H,y^Q(M; H) and 

oo 

\B * u\h^^^^^-h) < J e~''^\\B{t)\\ dt 



This completes the proof. □ 

Remark 3.2. Note that since B* commutes with do^i, we can extend B* to a bounded linear 
operator on H^^k{M; H) for each k £ 

Corollary 3.3. Let B G Li^^(IR>o; L(-ff)). Then limjy_j.oo \B\li ^(R-:^(y,L{H)) = o.iT'd thus, in 
particular liuii,^^ 

Proof. This is an immediate consequence of Lemma 13.11 and the theorem of monotone conver- 
gence. □ 



^Note that scalar analyticity on a norming set and local boundedness is equivalent to analyticity (see |14l 

Theorem 3.10.1]). 
^The integral is defined in the weak sense. 
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3.1 Hyperholic-type equations 



Lemma 3.4. Let B £ Li^^(M; for some fi > and u £ H^fi{R;H) for v > fi. Then 

{Cu{B * u)) (t) = V2^B{t - ii/) iCuu{t)) 

for almost every t G M. 

Proof. The proof is a classical computation using Fubini's Theorem and we omit it. □ 

Prom now on let Hq, Hi be complex Hilbert spaces, A : D(A) C Hq — )• Hi be a densely defined 
closed linear operator, Bi, B3,Ci,C3 G Li^^(]R>o; L(i/o)) and ,62, -B4, C2, C4 G Li^^(IR>o; 
for some /i > 0. We consider linear material laws of the form 



M{z) 



+ z 



Qi (^ai(-iz-i)) (5i(-iz-i)) 

Q3 (C3i-iz-^)Y' PslM-iz-^ 



Q4 fC4(-iz-i)) ^pJB^{-iz-^] 



for z G -Be (^^) ^ j 1 where Pi, Qi are complex affine linear functions, such that Qi \^Ci (— i-z ^) 

is boundedly invertible for z G i?c 5^) and 2; 1— )• Qi{Ci{—iz~^))~^ Pi{Bi{—iz~^)) is bounded 
for every i G {1, ... 4}. Under these assumptions M defines a linear material law and we con- 
sider evolutionary problems of the form 

(^9o,.M [d^l) + ( f ) ) C/ = F. (7) 

According to Theorem 12.91 it suffices to check the solvability condition ([6]), which will be done 
in concrete cases in the forthcoming subsections. Throughout we set H := Hq © Hi. 



3.1 Hyperbolic-type equations 

In this subsection we consider the case -Psj-Ri = 0, Pi{x) = 1 + \/27rx, Qi,P2 = 1 and 
Q2{x) = 1 — V^vrx. Thus the material law reads as follows 



/1 + V2^C(-iz-i) \ 



for some C G Li^fj^(M>Q; L(Hq)), B G Li^^(]R>o; i^(i^i)), where // > 0. According to Corollary 
13.31 there exists /xq > M) such that \\B * \\l{H„ a{^;Hi)) < 1 for each > fiQ. To ensure that the 
function M defines a linear material law, we choose r := Note that the system ([7]) with 

M given by ([8]) and F = Q G i?,.,o(M; H) reads as 

K'ra-irO^(-v;))(:)KO' 



11 



3 Integro-Differential Equations 



In the special case that g = we obtain 

do{l - B*)-'^q = Av 

or, equivalently, 

q = Oq^I - B*)Av. 
If we plug this representation of q into the first line of Equation ([9]) we get 

9o(l + C*)v + A*dQ^{l - B*)Av = f 

which gives, by defining u := Oq^v 

9^(1 + C*)u + A*{1- B*)Au = f. 

A semi-linear version of this equation was treated in [3] for scalar- valued kernels, where criteria 
for the well-posedness and the exponential stability were given. Also in |22| this type of 
equation was treated for scalar-valued kernels and besides well-posedness, the polynomial 
stability was addressed. In both works the well-posedness (the existence and uniqueness of 
mild solutions) was shown under certain conditions on the kernel by techniques developed for 
evolutionary integral equations (see |21j). We will show that the assumptions on the kernels 
made in both articles can be weakened such that the well-posedness of the problem can still 
be shown, even for operator-valued kernels. 

To ensure the well-posedness of ([9]) , we have to guarantee that there exist ri , c > with ri < r 
such that for all z G Bc{ri, ri) : 

d\t z~^il - V2^B{-iz^^)y^ > c (10) 

and 

me z-^{l + V2^C{-iz~^)) > c. (11) 

Remark 3.5. One standard assumption for scalar- valued kernels is absolute continuity. In our 
case this means that there exists a function G £ Li_^(M>o; L{Hi)) for some > such that 



B{t) = j G{s) ds + B{0) (tGM>o) 



for the kernel B. 

For simplicity let us assume C = 0. Note that due to the absolute continuity, B is an element of 
Li^;y(R>o; L{Hi)) for each i/ > fi and we choose i/ large enough, such that \B\i^ ^(R>o;L(_H'i)) < 1- 
In this case (|lUp can be easily verified. Using the Neumann-series we obtain 



oo 

-\l - V2^B{-iz-^))-^ = z-^ [V2^B{-iz-^] 

k=0 

oo 



z"^ + z"^\/2^S(-iz"^) (V2^B{-iz^^) 

k=0 
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3.1 Hyperbolic-type equations 



for z G Be 2^;) • The Fourier-transform of B can be computed by 

B{-iz-') = z (d{-iz-') + ^5(0)) 

and hence we can estimate 

mez-^il - V2^B{-iz-^))-^ 

oo 

= D\e z^^ + fHe (^V2^G{-iz^^) + -B(O)) ^ (^V2^B{-iz-^] 

k=0 



> V 



> V 



l-|^/2^5(-i(.)-^)|^^(5^(_L,_L);i(^,)) 
l^lii..(ffi>o;i(^^i)) + 11^(0)11^(^/1) 



l'^lil,^(K>Oi-t(«"l)) + ll'^(0)lli'{«"l) 

-, 2^; . omce 

yields the assertion. 



for every z G (21;, 2^) • Smce ^(g^^,^(„^„ ^ II-S(0)||l(//i) as ^ 00, this 



In the case, when C and i? are not assumed to be differentiable in a suitable sense, the 
conditions (jlOp and pT]) are hard to verify. We now state some hypotheses for B and C and 
show in the remaining part of this subsection, that these conditions imply ()10p and (jlip . 



Hypotheses. Let T G Li^^(]R>o; where G is an arbitrary Hilbert space and /i > 0. 

Then T satisfies the hypotheses ([i]),® and (|IiH) respectively, if 

(i) for all t G M>o the operator T(t) is selfadjoint, 

(ii) for all s, t G M>o the operators T{t) and T(s) commute, 

(iii) there exists I'o > n such that for all t G M 

t3mf{t - iuo) < 0. 

Remark 3.6. 

(a) If T satisfies the hypothesis ([1]), then 

00 

Jmf{t - iuo) = I sm{-ts)e~'""T{s) ds = -3mf{-t - m) (t G M) 

v2vr 7 


and thus (jm]) holds if and only if 

Jmf(t-ii/o) < (tGM>o). 



(b) Note that in |22j and [H the kernel is assumed to be real-valued. Thus, @ and du]) are 
trivially satisfied. In |22| we find the assumption, that the kernel should be non-increasing 
and non- negative, i.e., T(s) > and T{t) — T{s) < for each t > s > 0. Note that these 
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3 Integro-Differential Equations 



assumptions imply 

{{e-^'Tit) - e-^'T{s)) x\x) = e-''\{T{t) - T{s)) x\x) + (e"^* - e-'''){T{s)x\x) < 
for every t>s>0,z^>0 and x £ G. Hence, we estimate for t > and x £ G 
{3mf{t - iuo)x\x) 

oo 

^ sm{-ts)e-''°'{T{s)x\x) ds 



/(2fc+l)f 2(fc+l)f 

[ sm{-ts)e~''"'{T{s)x\x) ds+ [ sm{-ts)e-''°'{T{s)x\x) ds 

\ 2fcf (2fc+l)y 
/ {2fc+l)f 

/ sin( -18)6-""' (T( s)u\u) ds + 

(2fc+l)f 



+ y sin(-ts-7r)e-'^°(^+T) (|r(s + ^ 



u 



u) ds 



2fc2 



{2fe+l)^ 



1 /" 



'T(s + j]]u 



u) ds 



2fcf 
{2fc+l)f 



u) ds < 0, 



which yields according to (a). The authors of j3] assume that the integrated kernel 
defines a positive definite convolution operator on L2(IK>o)- However, according to [H 
Proposition 2.2 (a)], this condition also implies (jm]) . 

Lemma 3.7. Assume that T € Li^^(M>o; i^(G')) satisfies the hypotheses ^ and /fm|]. T/ien 
uie /lawe /or a// > 1^0 '^'^t^ t G M 

t3mf{t - ii/) < 0. 



Proof. Let x £ G and v > fi. We define the function 

/(t) := (T(t)x|x) (tGM) 
which is real-valued, due to the selfadjointness of T{t) and we estimate 

|/(t)|e-'^*dt< / ||r(t)||e-^*dt|x|2 
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3.1 Hyperbolic-type equations 



which shows / G Li^^(M>o). We observe that 



T{t — w)x\x 



e-'^'e-^'Tis) dsx 



1 



'2^ J \ V 7 I / 

= f{-t-iv) 

for each t G M, > /i. Hence, by 

{3mf{t - \v)x\x) = '2m{f{t - iv)x\x) = 3mf{-t - ii/) 

it suffices to prove 1 3m f{—t — ii/) < for v > vq, t G R under the condition that 1 3m f{—t — 
i'^o) ^ for each t G M. For this purpose we follow the strategy in [Sj Lemma 3.4] and employ 
the Poisson formula for the half plain (see |24l p. 149]) in order to compute the values of the 
harmonic function 3m /(•) : [M] — i [M>^] — )• R. This gives, using 3m /{—s — iu) = — 3m f{s — iv) 



1 



3m f{—t — \v) = — 



J {t-s)^ + {iy- i/o) 



2 3m /(— s — ii^o) ds 



TT 17 (t + s)2 + (i/ - i/o)2 



3mf{s — wq) ds + 



+ 



(t - sr + (z^ - i^o)'- 



3m f{—s — ii^o) ds 



1^-1^0 



1 



4t 



^ J V(t-s)2 + (l/-Z.o)2 (t + S)2 + (z. - I.o)2 


oo 

^' — T^Q f f S 



vr J y((t_5)2 + (^_^Q)2)((i + ^)2 + (^_^Q)2) 




3m f{—s — iuo) ds 
3m /(— s — iuf)) ds, 



which implies 

t3mf{—t — iu) 



= 4t 
< 0. 



2 1^-1^0 



J \{{t- S)2 + 1^0?) {{t + s)2 + (l. - 1.0)2) 





Jm/(— s — ii^o) ds 



□ 



Lemma 3.8. Let B satisfy the hypotheses ^-lin^. Then there exists < ri < r such that for 
all z G Bc(ri,ri) the condition U0\) is satisfied. 
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3 Integro-Differential Equations 



Proof. Let x £ Hi and set ri := min | • Let z £ B£{ri,ri) and note that z ^ = it + 

for some t G M.,u > i/q. Since the operator 1 — \/2tt B{t — \v) is bounded and boundedly 
invertible, so is its adjoint, which is given by 1 — \/^B{—t — iu) since B{s) is selfadjoint for 
each s G M. We compute 

$He((it+z^)(l-\/2^S(t-ii/))-^x|x) = d\t{{it+i^){\l-V2^B{t-w)\'^)-^{l-V2^B{-t-iu))x\x). 

We define the operator C := |1 — y/27rB{t — Furthermore, note that due to the 

assumption that the operators pairwise commute, we have that the operators B{-) commute, 
too. This especially implies, that B(t — iu) is normal and hence C and 1 — \/^B{—t — iu) 
commute. Thus, we can estimate the real part by 



9^e((it + i^)C'^{l - V2^B{-t - w))x\x) 
-- me{it + - V2^B{-t - iu))Cx\Cx) 

- V2^mcB{-t-iu)^ Cx\Cx) +t{V2^ 3m B{-t-iu)Cx\C x) 

^B{-t - \v)\\) \Cx\^ - V2^{t3m B{t - w)Cx\Cx) 



> 1/(1 - \B\L^^^^f,^^.L{Hi)))\Cx\'^, 
where we have used Lemma [3.71 Using now the inequality 

|x| = \C~^Cx\ = |(l - V2^B{t - ii/)) Cx\ < (1 + |SU,,,(R>o;L(Hi)))l^^^l 
we arrive at 

1 - |-B|li,,(R>o;L(//i)) I ,2 



D\t{{it + z^)C^(l - V2TTB{-t - ii/))x\x) > u- 



1 + I^Ili,,{R>o;L{Hi)) 



1 + l-^lLi,.o(iR>o;i(^^i)) 



which shows the assertion. 



□ 



After these preparations we can state our main theorem. 



Theorem 3.9 (Solution theory for hyperbolic-type integro-differential equations). Let A : 

D{A) C Hq —7- Hi he a densely defined closed linear operator and C G Li^^(M>o; I/(i?o)), -B G 
Li^^(M>o; L(i7i)) for some /i > 0. Assume that C and B are absolutely continuous or that B 
satisfies the hypotheses li))-/ fm)] and C satisfies ^ and / fm|] . Then the problem ^ is well-posed 
in the sense of Theorem \2. 91 



Proof. In case of absolute continuity of i? or C the positive definiteness condition is satisfied 
according to Remark 13.51 If B satisfies ()H)- ()iiI|) then Lemma 13.81 yields (jlOp . If C satisfies (P 
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3.1 Hyperbolic-type equations 



and dm]) then we estimate 

(it + I/) (l + V2^C{t - w)) =v{l + \/2^9^e C{t - ii/)) - \/2^t Jm C{t - iv) 



>v{l- |C|ij^(IR^(,.L(H))) 



for each t S M and v > uq, where we have used Lemma 13.71 Using Coronary 13.31 this yields 
the assertion. □ 

In apphcations it turns out that Equation ([9]) is just assumed to hold for positive times, i.e. 
on M>o and the equation is completed by initial conditions. So for instance, we can require 
that the unknowns v and q are supported on the positive real line and attain some given initial 
values at time 0. Then we arrive at a usual initial value problem. Since, due to the convolution 
with B and C the history of v and q has an influence on the equation for positive times, we 
can, instead of requiring an initial value at 0, prescribe the values of v and q on the whole 
negative real-line. This is a standard problem in delay-equations and it is usually treated by 
introducing so-called history-spaces (see e.g. |13| |9]). However, following the idea of |15| we 
can treat this kind of equations as a problem of the form ([9]) with a modified right-hand side. 
Let us treat the case of classical initial value problems first. 

Remark 3.10 (Initial value problem). For {f,g) G H^fi(M; H) with supp/, supp<^ C ]R>o we 
consider the differential equation 

*('ra-V.)^(-V„-))(:)^(^ 

on R>o completed by initial conditions of the form 

v{0+) = G D{A) and q{0+) = q^°^ G D{A*). 



We assume that the solvability conditions (jlOp and (jlip are fulfilled. Assume that a pair 



ft 

{v,q) G XR>o("^o)[-f^!^,i(I^) -f^)]!! solves this problem. Then we get 

. (1 + C* \fy-^^^^^vm\ / A*\(v\_(f 

on M, which is equivalent to 

(12) 

We claim that this equation is the proper formulation of the initial value problem in our frame- 
work. According to Corollarv 12 . 1 01 this equation admits a unique solution (v, q) G //j/^_i(M; H) 
and due to the causality of the solution operator we get suppv, supp q C M>o. We derive from 
(fT2]) that 



' (1 - B*r^ ) + [-A0 ))[q- XM>o ^ J " U J " ''^^^ V -Av^'^ 



^This means that we find a pair {w,p) £ H,^.i(R; H) that coincides with {v,q) for positive times. 
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3 Integro-Differential Equations 



which gives ( ^ ^^>o ® ^^^^ j g H,y^Q{M.; H). However, we also get thal| 



XK>o ® Q 

7\ f A* \ fv 



1 + c* Wi^-»>o®^^°^ 

(l-i?*)-V U-XR>o®9^°^ 



and hence, 



(0) ' 

v-XR^o'^v^^^ ) Gi/.,i(M;i?„i(|A*|+i)ei7_i(|^|+i)). 



!;>o ® 9' 



Using the Sobolev-embedding Theorem (see |20| Lemma 3.1.59] or |15l Lemma 5.2]) we obtain 

f V — Y]R (B) v^^^ \ 

that I /n> is continuous with values in H-^ i\A* \ + i) © H-\ (\A\ + i) and hence 

Q=(y- XR>o ® (0-) = (v- XR>o ® ^^^°^) (0+) 
in H_i{\A* \ + i) and thus 

^;(0+) =t>W in +i). 

Analogously we get 

g(0+) = g(0) in/?_i(|^|+i). 



Remark 3.11 (Problems with prescribed history). For (/, g) G iJiy^o(I^j H) with supp /, supp g C 
M>0 we again consider the equation 

*(^ra-V.)-(-v;))(:)^(0 

on R>o and the initial conditions 

7;|r<o = ^^(-oo), ^^(0+) = t;(_oo)(0-) and q\^^^ = g(_oo), ^(0+) = g(-oo)(0-). 

We assume that ^^(-oo) ^ Hi^^q{M.; Hq) with suppi;(_oo) C M<o and such that t'(„oo)(0— ) G 
D(^) and (l + C*)u(_oo) G XR>o("^o)[-f^v,i(^^; -^^o)] as weh as g(_oo) G i/,.,o(IK; -H'l) with 
suppg(_oo) C M<o, g(-oo)(0-) G D{A*) and (1 - S*)"^g(_oo) G XR>o(?"o)[^<.,i(M; ^^i)]- We 
want to determine an evolutionary equation for w := XR>o'^ P '■— XR>o9- We have that 

f\ / fi + c* V M ^* 

J - ^^>"PV {l-B*)-')^[-A 

°v ii-B*)'^) + [-A u+9{-oo) 



1 + c* \ ( A*\\ fw 
(l-i3*)-i ] + i j j U 



'Note that A and A*can be extended to bounded operators j4 : Ho —>■ H-i{\A*\+i) and A* : Hi — > _ff_i(|A|+i) 
respectively (cf. |20l Lemma 2.1.16]). 



18 



3.2 Parabolic-type Equations 



Hence, we arrive at the following equation for {w,p) : 



(14) 



. /1 + C* \,( A*\\(w\_(f\ ( il + C*)v^^^) 

{l-B.r^) + [-A ))[p)-[9)-^^>''^'[il-B*y\_^, 

Note that we can omit the cut-off function on the left hand side due to the causality of the 
operators. The conditions v{0+) = 'y(_oo)(0— ) and q{0+) = 9(-oo)(0~) now classical 
initial conditions for the unknowns w and p. Hence, following Remark 13.101 we end up with 
the following evolutionary equation for {w,p): 

l + C* \ /^(_)(0-) 

(i-i?*ri;^^U(-oo)(o-) 

(15) 



+ 



This equation possesses a unique solution in Hiy-i{M;H) with supp suppp C M^g due to 
the causality of the solution operator. Like in Remark 13.101 we get that 

--^«>o7(--)(0-))ei/.,o(M;i/) 
P-XR>o «'9(-oo)(0-) J 



from which we derive, using Equation ()15p . that 

w{0+) = ?;(_oo)(0-) and p{0+) = g(_oo)(0-) 

in i?_i(|A| +i) and //_i(|j4*| +i) respectively. We are now able to define the original solution 
by setting 

V ■.= w + ^(-oo) G H^fi{R; Ho) and q = p + g(„oo) G H^^o{R; Hi). 

Indeed, v and q satisfy the initial conditions by definition and on M>o the solution (v, q) 
satisfies the differential equation (fT3|) according to the computation done in ([K 



3.2 Parabolic-type Equations 

In this subsection we treat the case Pi{x) = 1 + \/27rx, Qi = l,P2 = 0,^3 = 0, P4 = 1 and 
Q4{x) = 1 — \/27rx. Hence, we end up with a linear material law of the form 

Miz)=(' + ^^^-''-'^')+z^^(\ " ,.-0 i^^Bcir,r)) 

\ oy [0 (1- V^B{-iz-^)^ ' ^ ^\ ' J) 



where B G Li^^(M>q; L{Hi)) and C G Li^fj_{M>Q; L{Hq)) for some /U > and r > is chosen 
suitably. The corresponding integro-differential equation is given by 



Ol \0 {1- B*)^' J \-A J J \q J \g 
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3 Integro-Differential Equations 



which, in the case g = 0, yields the parabohc equation 

dou + C* dou + A*Au - A* {B * Au) = f. (17) 

This problem was already considered in [3] for scalar- valued kernels. As it turns out the 
solution theory for this kind of problem is quite easy in comparison to the solution theory for 
the hyperbolic case. 

Theorem 3.12 (Solution theory for parabolic-type integro-differential equations). Let A : 

D{A) C Hq —7- Hi be a densely defined closed linear operator and B G Li^^{M.>q;L{Hi)),C G 
Li^^{M; L{Hq)) for some ^ > 0. Assume that C satisfies the hypotheses (0) and / fm|) (see 
Subsection 3.1). Then problem ligj) is well-posed in the sense of Theorem \2.9i 

Proof. To verify the solvability condition ([6]), we have to show that the operators y{ez~^{l + 
y/27^C{—iz~^)) and lHe(l — ^/27rB{—iz~^))^^ are uniformly strictly positive definite on some 
ball Bc{r,r). The first term can be estimated as in Theorem 13.91 For the second term we use 
the Neumann-series and estimate 

— 1 °° k 

mt (l - V2^B{-iz-^)^ =l + mc V2^B{-iz-^) (y2^B{-iz-^)^ 

k=0 

^ ^ _ SUp^eBc(r,r) II ^^^-^(-i^'^) IIl(Hi) 

1 - sup^gB^(,.^^) \\V2^B{-iz~'^)\\L{H^) 
^1 (r^0+), 

which yields the strict positive definiteness of the second term. □ 

Remark 3.13. In |27| the following kind of a parabolic-type integro-differential equation was 
considered: 

dou + C*dou + Au-B*Au = f, (18) 

where A : D{A) C Hq — t- Hq is a selfadjoint strictly positive definite operator and B,C G 
Li^^i(M.>q; L{Hq)). Equivalently we can consider the equation 

dou + (1 + C7*)-i(l - B*)Au = (1 + C*y^f. 

In this case the well-posedness can be shown, without imposing additional hypotheses on C. 
First we write the problem in the form given in ([7]). For doing so, consider the operator 
A : D{A) C Hi{^/A) — )• Hq. We compute the adjoint of this operator. First observe that for 
g G Hi{\^) we get 

{g\Af)Ho = {VAg\VAf)Ho 

for each / G D{A) and thus g G D{A*). Furthermore, if 5 G D{A*) there exists h G Hi{\^A) 
such that for all / G D{A) 

{g\Af)Ho = {h\f)H,i^^) = {^/Ah\y/Af)H, = {h\Af)H,. 
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Since A has dense range we conclude that g = h G Hi(^/A). Thus, the adjoint is given by the 
identity 1 : Hi{\/^) C Hq — )■ Hi{^/A). We rewrite Equation (|18p in the fohowing way 

^ ^\\(^\ = ({i + c*y'f\ 

{I- B*y\l + C*) J ^ \-A0 J J \q J \ J- 

Note that this is now an equation in the space Hiyfi{M.; Hi{y/A) © Hq). The strict positive 
definiteness of 9\e{l — B*)^^{l-irC*) follows from the strict positive definiteness of 9^e(l— 
(compare the proof of Theorem I3.12p and the fact that ||(1 — B*)~^C * ||l(_H'„ o(lK;^^o)) ~^ ^.s 
— )• oo. 




4 Application to Visco-Elasticity 

In this section we apply our findings of Subsection 3.1 to the equations of visco-elasticity. For 
doing so we first need to introduce the involved differential operators. Throughout let C M"- 
be an arbitrary domain. 

Definition 4.1. We consider the space 

L^inr^"^ := = {^^J)^J^{^_n} I Vi, j G {1, . . . ,n} : G L2(^^)} 
equipped with the inner product 

(^-1$) := J trace(^'(x)*$(x)) dx (^,$ G L2(0)'^^'^). 
a 

It is obvious that ^2(0)"^" becomes a Hilbert space and that 

H,y^{n) := £ L2((7)"^" I ^{xf = ^{x) (xGn a.e.)} 

defines a closed subspace of L2(0)"'^'^ and therefore Hsym{^) is also a Hilbert space. We 
introduce the operator 



Grad \c^^n)^ : C^in)" C L2(17)" ^ H,y^{n) 

{^i)ie{i,...,n} ^ (^^{9i(pj + dj^i) 



i,je{l,...,n} 



which turns out to be closable and we denote its closure by Grade . Moreover we define 



Div := - Grad* 



For $ G Cc(ri)"^" n i7sym(^^) one can compute Div$ by 



n 



I «£{!,. ..,n} 
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4 Application to Visco-Elasticity 



The equations of linear elasticity in a domain 0, C read as follows (see e.g. |10| p. 102 ff.]) 

do{gdou)-DwT = f (19) 
r = CGradcn, (20) 

where u € Hiyfi{M; L2{^)^) denotes the displacement field and T £ Hi,^q{M; Hsy^{Q)) denotes 
the stress tensor. Note that due to the domain of the operator Grade we have assumed an 
implicit boundary condition, which can be written as 



u 



on dn 



in case of a smooth boundaryjf] The function g G Loo(J^) describes the pressure and is assumed 
to be real-valued and strictly positive. The operator C G L(ffsym(^))> linking the stress and 
the strain tensor Grade it is assumed to be selfadjoint and strictly positive definite. In viscous 
media it turns out that the stress T does not only depend on the present state of the strain 
tensor, but also on its past. One way to model this relation is to add a convolution term in 
(I20D, i.e., 

t 

T{t) = CGiadcu{t)- j B{t-s)Giadcuis) ds, (t G M), (21) 

— oo 

where B G Li^^(IR>o; i^(-ffsym(f^)))- If we plug (|2ip into ()19p we end up with the equation, 
which was considered in [B] under the assumption, that B is absolutely continuous. 
We now show that (jlOp and ()2ip can be written as a system of the form Q . For this purpose we 
define v := Bqu. Note that the operator C — B* = C2 (l — C~2 (B*) C~2 J C2 is boundedly 



invertible, since C is boundedly invertible and yt C 2B{t)C ^ j ^ ^i,^t(IK>o; -^(^sym(^))) 

which gives that ^1 — C~2 (B*) C~2^ is boundedly invertible on H,^^q{W; ffsym(^)) for large 
I' (Corollary 13. 3p . Therefore we can write Equation (|2ip as 



and by differentiating the last equality we obtain 



(B*) C'^] ^ C-It = C^ Graded. 



Thus, we get formally 



1 



ofl-C-^i?*) ^^-^r'j^Ut^^ Grade ' " ^ j ) ( C^T ) " ( ' ^ 




(22) 

Following the strategy used in [20^ Section 5.2.1] we define a new inner-product on H in the 
following way. 



^By suitable realizations of the operators Grad and Div one can model also more complicated boundary 
conditions (see [251 p. 98 fT.]). 
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Definition 4.2. Let H he a. Hilbert space and M G L{H) be selfadjoint and strictly positive 
definite. Then we define an inner product on H b}0 

{Ay)M[H\ ■= {M~^x\y)H, {x,y G H) 
which yields an equivalent norm on H. We denote the Hilbert space (^H, {■\-)m[H]) by M[H]. 
Lemma 4.3. We define the operator 

A := Grade : i:>(Gradc) C g^^ [L2(J^)"] ^ Hsy^{n). 
Then A* = -^--^DivCi 

Proof. Let ^ £ D ^DivCz^ . Then we compute for ip £ D(Gradc) 

= (-01 - Q^'^B[vC^<P)^-i[L,(n)n], 
showing that -g^-^ DivC^ C A*. If $ G D{A*), then there exists t] £ ^2(5^)" such that 

The latter yields that £ L'(Div) and — Div Ca $ = gr], which gives the remaining operator 
inclusion. □ 

The previous lemma shows that (|22p as an equation in H^^q (M; ^-i[L2(f^)"] i?sym(f^)) is of 
the form given in ([9]). Thus our solution theory (Theorem 13. 9p applies. So, in the case that 
B is absolutely continuous (as it was assumed in |8]) we get the well-posedness. If we do not 
assume any smoothness for B we end up with the following result. 

Theorem 4.4. Assume that B satisfies the hypotheses B\)- [m\) and that C and B{t) commute 
for each t G R. Then \2'J3(l is well-posed as an equation in H^^ (M; ^?~^[L2(il)"] © //sym(f^)) 
for V large enough. 



Proof. Note that C 2 (i?*) C 2 = | C '^B{-)C 2 J *^ so according to Theorem 13.91 it suffices 



to verify the hypotheses (p|- (liiH) for the kernel C 2B(-)C 2. The conditions (ji]) and (jn]) are 
obvious. Furthermore 



C-2B{-)C-2 j {t - iu) = C~2B{t - \v)C—2 
for each i G M, > /i and hence by the selfadjointness of C 

Jm ( C^5B(.)C^I ) (t - ii/) = C7-i '2xnB{t - \v)C 



''Note that A/ ^ £ I^{H) is also selfadjoint and strictly positive definite. 
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= (tJm5(t-iz^o)C"2$ 

-ffsym(fi) 



1 



C~2$) <o 



This gives 

t3m (^C~^B{-)C~^^ {t - ii/o)$ 
for each <I> G Hsym{^) and t G M. □ 

5 Conclusion 

We have shown that hnear integro-differential equations with operator- valued kernels of hy- 
perbolic and parabolic type are covered by the class of evolutionary equations and we gave 
sufficient conditions for the well-posedness in both cases. Moreover, using the causality of 
the solution operators and its continuous extensions to the Sobolev-chain associated to the 
time-derivative, we have proposed a way to treat initial value problems without introducing 
history spaces (see Remark 13.101 and Remark I3.1ip . Note that most of the results can also 
be formulated for differential inclusions, where we replace the operator A* A by a maximal 
monotone relation, using the techniques developed in |26| I25|. 

So far, we have restricted ourselves to the case of affine linear functions Pi,Qi in ([8]). However, 
one could also treat polynomials instead, which would yield integro-differential equations with 
compositions of convolution operators. Also one could allow non-vanishing off-diagonal entries 
in the block operator matrices in ([8|). For these cases it is left to state some sufficient condi- 
tions for the well-posedness, which are easier to verify than the abstract solvability condition 

dSD. 

Another possible generalization would be integro-differential equations with convolutions with 
unbounded operators, for example equations of the form 

{di + A*{l-B*)A)u = f, 

where the operators B{t) are unbounded, but continuous as operators on the Sobolev-chain 
associated with the operator \A*\ + i. Also in this case the solution theory for evolutionary 
equations is applicable, if we assume that the resolvent (1 — B*)^^ is bounded. However, the 
main problem remains to give sufficient conditions for the solvability condition ([6]) in this case. 
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